Measurement of the Nucleon Structure Function F2 in the 
Nuclear Medium and Evaluation of its Moments 



M. Osipenko'J''^f'*, G. Ricco^, S. Simula^", M. Ripani^, M. Taiuti^, K.P. Adhikari'^^, 
M.J. Amaryan'*'', M. Anghinolfi'^ , H. Avakian'^^'P, H. Baghdasaryan'^\ M. Battaglieri'i, 
V. Batourine'^^, I. Bedlinskiy', A.S. Biselli'"'^'", D. BranfordJ, W.J. Briscoe", 
W.K. Brooks^j'^i^, V.D. Burkert^'^, S.L. Careccia^^, D.S. Carman^'^, P.L. Cole°'^'^, 
P. Collins^'i, V. Crede"", A. D'Angelo''^^, A. DanieP, N. Dashyan^", R. De Vita^, 
E. De SanctisP, A. Deur^*^, B Dey°, S. Dhamija', R. Dickson^, C. Djalali^s^ 
^ ■ D. Doughty'^''''\ R. Dupre'', H. Egiyan'^'^'", A. El Alaoui'', P. Eugenio'", S. Fegan'^'^, 

C\l ; T.A. Forest°'^^, A. FradP, M.Y. Gabrielyan\ N. Gevorgyan'^" , G.P. Gilfoyle^^ 

K.L. Giovanetti", W. Gohn', R.W. Gothe'^s^ k.A. Griffioen™, L. Guo'^*''^, K. Hafidi'^, 
5J; H. Hakobyan^j-'^", C. Hanretty™, N. HassalP'", D. Heddle''^'^'^, K. Hicks^ M. Holtrop'^, 

C/^ '■ Y. Ilieva^s^ d.G. Ireland^i^, E.L. Isupov'^f, S.S. Jawalkar™, H.S. Jo", K. Joo'-'^j, 

D. Keller^ M. Khandaker^, P. KhetarpaP^, W. Kim^, A. Klein^^, F.J. Klein^'^^^^, 
Psj ■ V. Kubarovsky'^'^ S.E. Kuhn'*'', S.V. Kuleshov^'j •* , V. Kuznetsov'^, K. Livingston'*'*, 

H.Y. Lu'^s^ D. Martinez", M. Mayer'*'*, J. McAndrewJ, M.E. McCracken", 

B. McKinnon'*'*, C.A. Meyer°, M. Mirazita?, V. Mokeev'*^''*^, B. Moreno^, K. Moriya^ 
(]J [ B. Morrison'*, H. Moutarde^, E. Munevar**, P. Nadel-Turonski^'^, R. Nasseripour'*^''', 

S. Niccolai^ I. Niculescu"*'**, A.I. Ostrovidov™, R. Paremuzyan'*'*, K. Park'*^'"''^^ 
O ; S. Park***, E. Pasyuk'*'^, S. Anefalos Pereira?, S. Pisano% O. Pogorelko*, 

^ ■ S. Pozdniakov*, J.W. Price'^, S. Procureur^, Y. Prok**''^, D. Protopopescu'*'*, 

B.A. Raue'''*'*, G. Rosner'*'*, P. Rossi?, F. Sabatie^^'*'*, M.S. Saini'**, J. Salamanca", 

C. Salgado'', P. Saracco'^, R.A. Schumacher'', H. Seraydaryan****, Y.G. Sharabian**'*-'*'*, 
<N ; E.S. Smith'*'*, D. Sober'', D. Sokhan^^^^ g_g^ Stepanyan^, S. Stepanyan'*'*-'***, P. Stoler'*'*, 

S. Strauch'*8, D.J. Tedeschi'*s, S. Tkachenko'*'*, M. Ungaro', B .Vernarsky^ 
^ ; M.F. Vineyard'*'-'*", E. Voutier*, D.P. Watts'*'*'^ D.P. Weygand'*'*, M.H. Wood'^, 

. A. Yegneswaran'*'*, J. Zhang'*'*, B. Zhao''^ 

cn 
(N 

^ , '^California State University, Dominguez Hills, Carson, CA 90747 

^3 . '^Canisius College, Buffalo, NY 

Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
f Catholic University of America, Washington, D.C. 20064 
^CEA, Centre de Saclay, Irfu/Service de Physique Nucleaire, 91191 Gif-sur-Yvette, France 
, '^Christopher Newport University, Newport News, Virginia 23606 

• ''University of Connecticut, Storrs, Connecticut 06269 

i ^Edinburgh University, Edinburgh EH9 3JZ, United Kingdom 

'^Fairfield University, Fairfield CT 06824 
^Florida International University, Miami, Florida 33199 
'^Florida State University, Tallahassee, Florida 32306 
^The Ceorge Washington University, Washington, DC 20052 

"Idaho State University, Pocatello, Idaho 83209 
f INFN, Laboratori Nazionali di Frascati, 00044 Frascati, Italy 
'^INFN, Sezione di Cenova, 16146 Cenova, Italy 
^INFN, Sezione di Roma Tor Vergata, 00133 Rome, Italy 
"Institut de Physique Nucleaire ORSAY, Orsay, France 
^Institute of Theoretical and Experimental Physics, Moscow, 117259, Russia 
James Madison University, Harrisonburg, Virginia 22807 



X 



'^Argonne National Laboratory, Argonne, Illinois 6044^ 
^Arizona State University, Tempe, Arizona 85287-1504 



1 



'" Kyungpook National University, Daegu 702-701, Republic of Korea 
■^LPSC, Universite Joseph Fourier, CNRS/IN2P3, INPG, Grenoble, France 
^ University of New Hampshire, Durham, New Hampshire 03824-3568 
^Norfolk State University, Norfolk, Virginia 23504 

'''Ohio University, Athens, Ohio 4 5701 
"■"■Old Dominion University, Norfolk, Virginia 23529 
"''Rensselaer Polyteclm-ic Institute, Troy, New York 12180-3590 
"■"^ University of Richmond, Richmond, Virginia 23173 
"'^Universita' di Roma Tor Vergata, 00133 Rome, Italy 
"■^INFN, Sezione di ROMA HI, 00146 Rome, Italy 
"■^ Skobeltsyn Nuclear Physics Institute, Skobeltsyn Nuclear Physics Institute, 119899 Moscow, Russia 
University of South Carol-ma, Columbia, South Carolina 29208 
"■^Thomas Jefferson National Accelerator Facility, Newport News, Virginia 23606 
'"■Union College, Schenectady, NY 12308 
"■^ Universidad Tecnica Federico Santa Maria, Casilla 110-V Valparaiso, Chile 
"■''University of Glasgow, Glasgow G12 8QQ, United Kingdom 
"'University of Virginia, Charlottesville, Virginia 22901 
"■^College of William and Mary, Williamsburg, Virginia 23187-8795 
"^Yerevan Physics Institute, 375036 Yerevan, Armenia 



Abstract 

We report on the measurement of inclusive electron scattering off a carbon target per- 
formed with CLAS at Jefferson Laboratory. A combination of three different beam ener- 
gies 1.161, 2.261 and 4.461 GeV allowed us to reach an invariant mass of the final-state 
hadronic system W « 2.4 GeV with four-momentum transfers ranging from 0.2 to 
5 (GeV/c)^. These data, together with previous measurements of the inclusive electron 
scattering off proton and deuteron, which cover a similar continuous two-dimensional 
region of and Bjorken variable x, permit the study of nuclear modifications of the 
nucleon structure. By using these, as well as other world data, we evaluated the F2 
structure function and its moments. Using an OPE-based twist expansion, we studied 
the Q^-evolution of the moments, obtaining a separation of the leading-twist and the 
total higher-twist terms. The carbon-to-deuteron ratio of the leading-twist contributions 
to the F2 moments exhibits the well known EMC effect, compatible with that discov- 
ered previously in a;-space. The total higher-twist term in the carbon nucleus appears, 
although with large systematic uncertainties, to be smaller with respect to the deuteron 
case for n < 7, suggesting partial parton deconfinement in nuclear matter. We speculate 
that the spatial extension of the nucleon is changed when it is immersed in the nuclear 
medium. 

Keywords: moments, nuclear modifications, nucleon structure, higher twists, QCD, 
OPE 
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1. Introduction 

Inclusive electron scattering off a nucleus, e{k)+A{P) e'{k')+X, may be expressed, 
under the one photon exchange approximation, in terms of the total absorption of a 
virtual photon by the nucleus: 7*(q = fc — k') + A{P) -> X. At sufficiently large 
squared four- momentum transfer, = —q'^, the coherent electron- nucleus scattering 
contribution is negligible, facilitating a description of the scattering in terms of individual 
nucleons in the nuclear matter. The inclusive electron- nucleon scattering cross section 
is related to the parton momentum distribution along the virtual photon (longitudinal) 
direction. This structure, however, depends on the environment in which the nucleon is 
embedded. In particular, the nuclear medium environment interferes with the nucleon 
internal structure. The experimental evidence for such nuclear modifications in the 
Deep Inelastic Scattering (DIS) regime was found by the European Muon Collaboration 
(EMC) ^. Other similar observations came from intermediate energy experiments that 
observed changes in the nucleon and nucleon excited state form- factors [yl- In order to 
study the transition of the nuclear modifications from the DIS regime (EMC effect) to the 
larger spacial scale region (the nucleon and its excited state form- factors), the structure 
function F2 is calculated from the cross section and used to perform a Mellin transform to 
moments M'~'^ . The Mellin transform allows to express the observable from Minkowski 
x-space (x- representation of F2{x)) in Euclidean space of moments (n-representation of 
M^^), where Lattice QCD method works. 

The Cornwall-Norton moments M^'^ of the nucleon (in the nucleus with mass A-Ia) 
structure functions are defined as: 

M^^iQ^)^ dxa;("-2)i^2(a;,Q2), n > 2, n even, (1) 

JQ 

where M is the nucleon mass, x — -~Aq^/2Pq is the Bjorken variable per nucleon and n is 
the order of the moment. Quantum Chromodynamics (QCD) predicts the Q^-evolution 
of these moments for each given n Moreover, the DIS component of these moments 
and the contribution of the nucleon and the nucleon excited state form-factors have 
different behaviors Therefore, the study of the dependence of the structure 
function moments is important for understanding the nucleon modifications in the nuclear 
medium. 

The identification of these nuclear modifications to the nucleon structure in inclusive 
electron scattering is complicated by the presence of other nuclear effects S]. The most 
important effect comes from the motion of the nucleon in the nucleus, known as Fermi 
motion. There are also other effects that cannot be reduced to the total absorption of 
the virtual photon by the nucleon. These are shadowing. Final State Interactions (FSIs) 
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Figure 1: Experimental data on the carbon structure function F2{x, Q^) used for the moment evaluation 
in the CLAS kinematic region. The points show the world data from Refs. fiol. Hi]. HI. [13. [li. [13. \\K 
[13, [H, [ii, 113, m, [13, IM in Hi, [23, Hi, m, [li, [13, m, m m HI . The shaded area shows the CLAS 
data region. 



and Meson Exchange Currents (MECs). However, the latter ones are typicaUy smaher 
and can be locahzed by their kinematic regions of influence. 

The study of nuclear modifications can be performed by comparing the nucleon struc- 
ture of a free nucleon to one bound in the nuclear medium. To this end, the proton and 
deuteron structure function moments were previously measured in Refs. [a, Hi- The 
present article maintains the same framework of these two previous measurements in 
order to avoid additional bias due to the analysis technique. 

Similar studies were performed in the past in Ref. [8j and recently repeated by the 
Hall C Collaboration at Jefferson Lab [§], and reported contradictory results. In this 
paper, we increase the precision of the previous studies by a new measurement of un- 
polarized inclusive electron scattering on carbon, performed with the CLAS detector in 
Hall B at Jefferson Lab. Like in the case of the proton and deuteron measurements, the 
new data cover a wide continuous two-dimensional region in x and as shown in Fig. [T] 

In Section[2]we describe the QCD picture of the F2 moments. In Section|3]we discuss 
some data analysis details. For further information we refer to Refs. 0,01, from which 
many procedures are adopted. Finally, in Sections [3| and [S] we present our main results. 

2. QCD Analysis of Nucleon Structure 

QCD-based descriptions of inclusive lepton scattering can require Operator Product 
Expansion (OPE) machinery (3| . The OPE expresses the structure function moments in 
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a power series: 



M^'^m = Enr{f^r,t^f,asn,rMi^^j , (2) 

where k — 1,2,...,cxd, /i/ (fir) is the factorization (renormahzation) scalc"'^ and /i is an 
arbitrary reference scale. OnrCM/) is the reduced matrix element of the local operators 
with definite spin n and twist t (dimension minus spin). This quantity is related to the 
partonic structure of the target. Enrif^r, l^f,cts) is a dimensionless coefficient function. 
For sufficiently large and therefore small running coupling constant as{Q^), it can be 
analytically calculated in perturbative QCD (pQCD). 
The first term in this expansion (the leading twist): 

VniQ^) = £^n2(Mr,M/ja5)0„2(A*/), (3) 

can be completely evaluated in QCD. The operator matrix elements On2, given by the 
expectation values of the parton momentum distributions in the nucleon, are calculated 
in Lattice QCD 37] for n = 2, 4. The coefficient functions are calculated in pQCD at 
a'g precision, including resummation of specific effects in the large-x region to In^ (o^s) 



precision [36|, |39|. This term is responsible for the DIS inclusive cross section and its 
modification in nuclei gives rise to the EMC effect. 

Further terms are far less known, but can be studied phenomenologically by exploiting 
the general form of their Q^-dependence. The separation of individual twists related to 
specific QCD operators is not possible without knowledge of the corresponding coefhcient 
functions, but the total higher-twist term: 

2 °° / 2 \ 5 ('''—4) 

HTniQ^) = ^J2Er.r{t^r,H,C.s)0.n,{iif)i^^j (4) 

contribution can still be analyzed. This total higher-twist term is responsible for the 
multi-parton correlations in the nucleon and, therefore, it is connected also to the nucleon 
elastic form-factors and to the transition form-factors to the nucleon excited states. If 
there are nucleon modifications of these properties that are important at lower Q^, higher 
twists can be used to investigate them. 

In summary, a comparison of the leading-twist moments of the nucleon structure 
function in complex nuclei to that in the deuteron allows a study of the EMC effect in 
the moment space (M„((5^)). This allows a direct comparison to Lattice QCD calcu- 
lations of the nucleon modified by the presence of the medium. A similar comparison 
between the total higher- twist moments in deuteron and nuclei allows study of the nuclear 
modifications of confinement and multi-parton correlations. 

3. Data Analysis 

The data were collected at Jefferson Lab in Hall B with the CEBAF Large Acceptance 
Spectrometer (CLAS) using a solid, 0.18 g/cm^ thick, carbon target during the electron 



^We are working in the Soft Gluon Re-summation (SGR) scheme |39| . where fij = = Q 



beam running period in April-May 1999. The average luminosity was 6x10^'^ cm~^s~^. 
Data were recorded using electron beam energies Eq = 1.161, 2.261 and 4.461 GeV 
to maximize the and x kinematics range. The accumulated statistics at the three 
energies allows a measurement of the inclusive cross section with a statistical uncertainty 
(< 5%) in X and intervals of Ax =0.009 and AQ^ =0.05-0.1 (GeV/c)2. 

CLAS is a magnetic spectrometer [i^ based on a six-coil torus magnet whose field is 
primarily oriented along the azimuthal direction. The sectors, located between the mag- 
net coils, are individually instrumented to form six independent magnetic spectrometers. 



The particle detection system includes Drift Chambers (DC) for track reconstruction [4l| , 
Scintillation Counters (SC) for time-of- flight measurements [42"|, Cherenkov Counters 
(CC) for electron identification [4^, and Electromagnetic Calorimeters (EC) to measure 
neutrals and to improve electron-pion separation [4J|. The EC detectors, which have a 
granularity defined by triangular cells in a plane perpendicular to the incoming particles, 
are used to study the shape of the electromagnetic shower and are longitudinally divided 
into two parts with the inner part acting as a pre-shower. 

Charged particles can be detected and identified for momenta down to 0.2 (GeV/c) 
and for polar angles between 8° and 142°. The CLAS superconducting coils limit the 
acceptance for charged hadrons from about 80% at 6* = 90° to about 50% at forward 
angles {9 = 20°). The total angular acceptance for electrons is about 1.5 sr. Electron 
momentum resolution is a function of the scattered electron angle and varies from 0.5% 
for 6 < 30° up to 1-2% for 6 > 30°. The angular resolution is approximately constant, 
approaching 1 mrad for polar and 4 mrad for azimuthal angles: the resolution for the 
momentum transfer ranges therefore from 0.2 up to 0.5 %. The scattered electron missing 
mass (W) resolution was estimated to be 2.5 MeV for a beam energy less than 3 GeV 
and about 7 MeV for larger energies. To study all possible multi-particle states, the 
acquisition trigger was configured to require at least one electron candidate in any of the 
sectors, where an electron candidate was defined as the coincidence of a signal in the EC 
and Cherenkov modules for any one of the sectors. 

The data analysis procedure has been described in detail in Refs. @, 01- Therefore, in 
this article we focus on changes and improvements in the analysis. The most important 
improvements, leading to a significant reduction of the estimated systematic uncertainties 
relative to those of Refs. i6|i7|], are described in the following sections. 

3.1. Generic Procedures 

Phenomenological corrections to the reconstructed charged particle momentum were 
applied to compensate for small drift chamber misalignments and torus magnetic field 
map inaccuracies. These corrections were determined by subtracting the missing mass 
spectrum (W) measured using a solid carbon target from a measurement using CH2 in 
order to extract the proton elastic scattering response. The scattered electron momen- 
tum measurement was shifted for the entire W spectrum in such a way that the elastic 
scattering peak was at the correct position oiW = 938 MeV/c^. This was accomplished 
by making a fit of the hydrogen elastic peak distributions for a few azimuthal angle bins 
inside each sector under the assumption that the polar angle was measured precisely. 
These corrections were applied for each sector separately and improved the elastic miss- 
ing mass resolution by 15%. For the highest beam energy data set, the statistics from 
hydrogen were insufficient to determine momentum corrections. However, the Fermi mo- 
tion at these values leads to a very smooth cross section having no peaks or other 
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Figure 2: Hydrogen elastic scattering cross section extracted by means of CH2-C subtraction for three 
different beam energies: full squares - 1.161 GeV, full triangles - 2.261 GeV, reversed full triangles - 
4.461 GeV. The inner error bars show statistical uncertainties, while the outer error bars represent the 
statistical and systematic uncertainties summed in quadrature. The cross sections are compared with 
the data from Ref. [Js^ (empty circles) and with parametrization based on dipole form-factors (dashed 
line) including also radiative corrections (solid line). 



structures and, therefore, makes it relatively insensitive to possible small systematic un- 
certainties in the measured momentum. To check the absolute normalization of our data 
we subtracted the carbon data from CH2 yield and obtained the cross section of the scat- 
tering on hydrogen. The elastic scattering cross section was selected and compared to the 
previous data from Ref. 145j and to the dipole parametrization. The comparison, shown 
in Fig. [5J exhibits fairly good agreement, although with large systematic uncertainties, 
related to the small fraction of hydrogen events in the total CH2 yield. 

Although the results of the present analysis use realistic Monte Carlo simulations of 
the CLAS detector, such simulations have limited reliability at the edges of the detector's 
sensitive volume. There are also small regions of the CLAS detector that were very 
inefficient and, in some cases, were malfunctioning. These regions were eliminated using 
kinematic cuts to simplify the analysis. The kinematic cuts were defined as functions of 
the measured scattered electron momentum and angle to establish fiducial volumes with 
a uniform detector efficiency in the polar angle 0. These cuts along the edges of the 
acceptance removed about 50% of the recorded data. 

3.2. Electron Identification 

Electrons were identified in CLAS by a coincidence between signals from the Electro- 
magnetic Calorimeter (EC) and the Cherenkov Counter (CC). Furthermore, only events 
with a reconstructed Drift Chamber (DC) negatively charged particle track that matches 
these hits were selected. 
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The CLAS Cherenkov Counter (CC) inefficiency does not exceed 2% within the fidu- 
cial regions of the CLAS acceptance [1, However, the electron- pion separation does 
not work at the same level of precision and some pion contamination shows up in the 
electron candidate sample. This is especially important at high beam energies (e.g. for 
the 4.461 GeV data set), where the number of low momentum pions is significant. 

The pion contamination in the electron candidate sample appears as a single photo- 
electron peak in the measured CC spectra (see Fig. [3]). The main source of this contam- 
ination is the result of an accidental coincidence between a pion produced by quasi-real 
photoproduction with random noise in one of the 36 CC PMTs. This coincidence takes 
place due to the loose matching between the CC hit and the track within one sector of 
CLAS. A procedure was developed in Ref. 0] to geometrically and temporally match a 
CC hit with a track seen in the DC and a hit in the SC. The resulting Cherenkov distribu- 
tion is shown in Fig. |31 The reconstruction algorithm eliminated the pion contamination 
in the single-photoelectron peak, which in the worst case could reach 20%. 

The Electromagnetic Calorimeter (EC) is used to separate pions from electron can- 
didates with momenta above the pion threshold (2.7 GeV/c). Electrons entering in the 
EC release 30% of their energy in the sensitive volume on average, while pion losses 
are constant (see Fig. |4]). This EC property was exploited at large particle momenta 
for electron-pion separation by selecting particles with an energy fraction released in the 
EC above 20%. More details on this procedure can be found in Ref. [6]. Furthermore, 
pions just above the Cherenkov threshold produce less Cherenkov light with respect to 
electrons of the same momenta and, therefore, can be removed by a cut on the number 
of photoelectrons measured in the CC. 




10 X Number of photoelectrons 

Figure 3: The spectrum of photoelectrons measured in one sector of the CC at i/ = 3.4 GeV (beam 
energy 4.461 GeV). The hatched area represents the CC spectrum after applying the matching procedure 
described in the text. 
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3.3. Background 

The creation of e+e" pairs by real photons or through hadronic decays and the 
scattering of electrons from materials other than the target are the dominant sources of 
background electrons. 

3.3.1. e^e^ pair production 

The main source of e'^e^ pairs entering the CLAS acceptance is from tt^ production, 
followed by either Dalitz decay to 76+6" or 7r° decay to 77, where one of the photons 
converts to an e+e~ pair. A model using a Wiser fit to inclusive pion photoproduction 
was used to effectively reproduce the positron cross sections observed in a previous CLAS 
experiment ^zj. The above model was used to estimate the e+e" contribution to the 
measured inclusive cross section. The value of the contribution was assumed to be 
equal to the ratio of the inclusive e+ production cross section over the inclusive nuclear 
scattering parametrization including radiative processes (tail from elastic peak, 

bremsstrahlung and Schwinger correction): 



^e+e-V-C'0,a;, y j — — g— ;7^^^ ^ — , (5j 



where cre+ is the inclusive e"*" production cross section, equal to the integral of the tt" 
quasi-real electroproduction cross section folded with its decay probability and inte- 
grated over the allowed kinematic range, and t*^^ is the inclusive nuclear scattering 
parametrization [1, H^l folded with radiative processes, as described in Ref. Q. 
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In this parametrization the tt^ quasi-real electroproduction cross section is given by: 



dae 



p2 



dphdVLh Eh Je 



dE^ (t. 



E-y 



Eh 



d^(j~, 



^hX 



(6) 



where cr-yp^hx is the parametrized photoproduction cross section, Ph and Eh are the 
pion momentum and energy, px is the pion transverse momentum with respect to the 
photon direction, Eq is the electron beam energy, E^ is the exchanged photon energy, 
tint and text are the internal and external radiator thicknesses and Ml = \/pt ~ ^ht 
with m/i being the pion mass. Ci and C2 are the Wiser fit parameters. 

The external radiator thickness is small compared to the second term, t^xt — 0.0053 
radiation lengths. The internal radiator length used to estimate the e"*" rate (previously 
assumed to be 5%) was evaluated according to the Weizsacker- Williams formula [iof : 



tint — ^ 
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log- 



Eo + E' 



(7) 



where E' is the scattered electron energy and nie is the electron mass. 

The largest e~^e~ contribution is less than 9% of the inclusive cross section at x of 
0.2 as shown in Fig. [5] and only impacts the data taken at 2.261 and 4.461 GeV beam 
energies, where such low values of x are accessible. 




0.35 



0.4 



0.45 0.5 



Figure 5: The contribution of e+e" pair production events in the inclusive cross section at Q'^ = 1.55 
(GeV/c)^. The numbers of positrons and electrons A^^^ were obtained from different runs with 

opposite CLAS torus magnetic fields, allowing for a cancellation of the acceptance and efficiency in the 
ratio. The points show the measured quantity 1 — A''^+ /N^- , which represents the ratio of the number 
of electrons inelastically scattered off carbon to the total number of measured electrons. The curve 
represents the calculations from Eq. [5] 
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3.3.2. Target Cell Contribution 

The target wall contribution to the measured cross section has been removed using 
empty target data in bins of the energy Eq, the Bjorken scaling variable x and such 
that: ^ 

Nc{Eo,x,Q^) = NfMEo,x,Q^)--^^N,^pty{Eo,x,Q^) , (8) 

where Nf^ii , Ngmpty and Nc are the number of events in the full target data, the empty 
target data and the number of events due to scattering off the carbon target, respectively. 
Q/^" and Qt^^*^ represent the charge accumulated in full and empty target runs as 
measured by a Faraday Cup (FC) gated by the data acquisition live time. Any inaccuracy 
of the Z-vertex reconstruction (see Fig. did not affect the extracted cross section 
because this information was not used in the evaluation procedure. 

Multiple scattering and elastic electron-nucleus scattering at the lowest beam energy 
of Eq ~ 1.161 GeV caused the incident electron beam to diverge after passing through a 
full target such that a small fraction of the beam was no longer within the FC acceptance. 
This resulted in charge loss. A precise experimental determination of the charge loss was 
performed using the measured cross section of electron scattering by the empty target 
walls (Mylar/ Aluminum windows of the empty cell for liquids). Electrons scattered by 
the empty cell walls were selected according to their reconstructed Z-vertex shown in 
Fig. [6] The evident enhancement of the full target data is due to the FC charge loss. 
The overall normalization factor was estimated by selecting the region of Z from -6 cm 
to 3 cm in both empty and full target data and evaluating the ratio: 



rz=3cm full , 
J z — — 6cm tot 
nz=3cm ^empty 
J z— — 6cm tot 



rp J z — — 6cm tot 

J'norm — ^=3cm empty . ' > 



where cr^""'^™P*2'(^) = ^ i A'^/"'itT"^'^(-^) is the total cross section for the target 

corresponding to a given Z interval (for definition of variables see Section |3.5|) . The 
correction factor was found to be E^-J^^Jf^'^^ — 1.1459 with negligible statistical uncer- 
tainty. Selecting only the first empty target cell window using the cut on Z from -6 
cm to -1 cm gives the same ratio within 0.2%. The corresponding Z-vertex distribution 
shows good agreement between the empty and full target cross sections using the above 
correction factor Enorm- Numerical estimates for Eq = 4.461 GeV data give the ratio 
Fnofrk^''^ = 1-0006 ± 0.006 for the Z-cut from -10 cm to -2.5 cm. A systematic uncer- 
tainty of 3% was assigned to the FC charge loss correction for the Eq — 1.161 GeV data 
set based on the above analysis. 

3.4- Acceptance and Efficiency Evaluation 

Determination of the acceptance and efficiency corrections was based entirely on the 
Monte Carlo (MC) simulations developed for CLAS. Moreover, the systematic uncer- 
tainties of these corrections were estimated from a comparison of MC simulations with 
experimental data using a realistic model in the event generator. In short, the proce- 
dure was the following: we generated events with the event generator describing elastic, 
quasi-elastic and inelastic eC-scattering processes including radiative corrections. These 
events then were processed with the GEANT-based CLAS software simulating the detec- 
tor response. After that the standard CLAS event reconstruction procedure was applied. 
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Figure 6: vertex distribution for the carbon target in (open histogram), and carbon target out (hatched 
histogram) at Eq = 1.161 GeV. The large peak that can be partially seen on the right side of the plot 
is due to scattering off the carbon plate, while the three smaller peaks are due to the empty target 
windows. 



Finally, the ratio of reconstructed events to the number of generated events gave a com- 
bined efficiency/acceptance correction in each kinematic bin. 

Detector response simulations were performed in the same way as described in Ref. Q . 
The following improvements and changes for electron-nucleus scattering were imple- 
mented: 

1. Electron scattering events were generated by a random event generator with the 
probability distributed according to cr*^^. The values for the elastic and inelastic 
cross sections for electron-carbon scattering were taken from existing fits of world 
data in Refs. [H^ and [s, 48|, respectively. 



The contribution from internal radiative processes was added according to calcu- 



lations 51 



3. The event rate obtained in the simulations was then compared to the data, pre- 
serving the original normalizations (accumulated FC charge for the data and the 
number of generated events over the integrated cross section of the event generator 
for simulations). These normalized yields differ from the cross section by the ac- 
ceptance, efficiency and radiative corrections. The simulated events were subjected 
to the same cuts applied to the data. e~^e~ pair production and empty target 
backgrounds were subtracted from the data. The normalized yields obtained with 
the same set of cuts from the data and simulations were compared and found to be 
in good agreement at a; < 1.5 (see Fig. [7]). At larger x, the contribution of nucleon 
short range correlations (SRC), which are not included in the model, may play a 
significant role 52, 53| . 

The efficiency appears to be fairly fiat and not lower than 97% inside the region of the 
detector defined by the fiducial cuts. 
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Figure 7: Comparison of normalized yields from the data (open circles) and reconstructed Monte Carlo 
(open triangles) for four different kinematic settings. The error bars include only statistical uncertainties. 



3.5. Structure Function F2{x,Q'^) 

The measured electron yield Nc, normalized to the integrated luminosity in conjunc- 
tion with Monte Carlo simulations, were used to extract the structure function in each 
kinematic bin. The Monte Carlo events were used to simultaneously obtain efficiency, 
acceptance, bin centering and radiative corrections. F2 has been determined using: 

p/.^^2^_ 1 J ^ Ncix,Q^) 2^ (-,r.^ 

t2\X,Ll ) — , i_e — 1 -7^^e+e-\X,Li ) , [W) 

p^ALQtot cTMott 1 + VlTfl *(^' ) 
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where v = Eq — E' \s the energy transfer, p is the density of the carbon target, Na is the 
Avogadro constant, Ma is the target molar mass, L is the target length, Qtot is the total 
charge in the Faraday Cup (FC) and ^'(a;, Q^) is the efficiency including the radiative 
and bin-centering correction factors: 

with: 



^ and vl,,„^iAr^_ (12) 



^'e/ / is the ratio between the number of reconstructed and generated events in the bin 
and a^^ is the event generator model cross section. The integral in Eq. [TH was taken 
over the bin area Ar. Here e is the virtual photon polarization parameter: 



tan^ - . (13) 



The Mott cross section g mou and the Jacobian J of the transformation between dVldE' 
to dxdQ^ are defined by: 



cos^ f , ^ xEf)E' 



4:E^ sm 



The structure function F2{x,Q^) was extracted using the fit of the function R{x,Q^) = 
'^l/o't described in Ref. [7\. However, the structure function F2 in the relevant kinematic 
range is very insensitive to the value of R. For example, at the typical kinematics of 
e — 0.75 and assuming a SLAC DIS R value of 0.18, the relative uncertainties of F2 and 
R are related by AF2/F2 = 0.03AR/ R. Therefore, at this kinematics a relative error in 
R of 20% will generate only a relative error of 0.6% in F2. The overlapped data from two 
different beam energies were combined using a weighted-average technique. Moreover, we 
checked that the parametrization used for R was consistent with the difference between 
the two cross sections within statistical and systematic uncertainties. 

Fig. [5] shows a comparison between the F2 data from CLAS and the other world data 
in a few bins. The CLAS data agree very well with all previous measurements. The 
measured F2{x,Q^) values are given in the center of the corresponding x and bin. 
The Monte Carlo based bin centering correction from Eq. [T^] allows to interpolate the 
data averaged over the bin to its value in the middle of the bin. Therefore, we do not 
provide any uncertainty in the determination of x and Q^. However, this uncertainty is 
implicitly included in the systematic uncertainties of Monte Carlo simulations. Indeed, 
Monte Carlo simulations include the effect of event migration from the original bin to 
neighbor bins, due to the finite detector resolution and energy loss in the target. Hence, 
the loss or gain of events in each single bin in the data is compensated by a similar change 
of the number of Monte Carlo events. As a test, we verified that the data obtained within 
significantly larger, not rectangular bins (in which bin migration is suppressed) agree with 
the presented data within statistical uncertainties. The values of F2{x, Of), together with 



their statistical and systematic uncertainties, are tabulated elsewhere 54 1. 
The radiative correction factor Vl/rad was calculated in the following way: 
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• the eC elastic radiative tail was calculated according to the "exact" Mo and Tsai 
formula (sH : 



• in the quasi-elastic peak region {W^^ + AW < W < 1.2 GeV) we have applied the 
correction formula to the continuum spectrum given in Ref. [5l| , which is based 
on the peaking approximation and is known to be reliable only when E' /Eq > 0.5. 
Here W'^'^ is the eC elastic peak position and AW is its width; 

• a.tW > 1.2 GeV we applied the exact Mo and Tsai formula to the quasi-elastic tail, 
while a formula based on the peaking approximation (referred to as the "unfolding 
procedure" ) was used for the inelastic spectrum. For an exact calculation of the 
quasi-elastic tail, it was necessary to extract quasi-elastic form-factors. To this end 
we integrated the quasi-elastic cross section given by our parametrization from the 
beginning of the peak up to = 1.2 GeV and performed a separation of the electric 
and magnetic form-factors. The tails calculated in these two kinematic regions do 
not exhibit any discontinuity at the point W ~ 1.2 GeV. This assured us that the 
peaking approximation formalism is safely applicable to the quasi-elastic tail up to 
W = 1.2 GeV. 

The radiative correction factor '^rad varies strongly in the explored kinematic range from 
0.7 up to 1.5. Fortunately, the largest corrections are given by the tails of the elastic and 
quasi-elastic peaks, for which calculations are very accurate (of the order of percent, see 



Refs. (5l|,l55l) 



3.6. Systematic Uncertainties 

Here we summarize all the systematic uncertainties taken into account in the present 
analysis: 

• The cross section overall normalization carries a systematic uncertainty due to the 
known precision of the target length and thickness. We estimated this uncertainty 
to be 1%. For the lowest beam energy set, Eq = 1.161 GeV, an additional system- 
atic uncertainty due to FC charge loss corrections of 3% was added. 

• Previous CLAS inclusive measurements 0| on proton and deuteron targets 
showed that the combined efficiency /acceptance systematic uncertainty of CLAS 
data averaged to 4.5%. We assumed a 3% uniform systematic uncertainty due to 
the efficiency evaluation, while the systematic uncertainty on the acceptance, based 
on Monte Carlo simulations, was evaluated separately. 

• There are two systematic uncertainties in the Monte Carlo simulation. The first 
one is due to the model dependence of the reaction cross section used for generating 
the events. Taking into account the good agreement between the data and Monte 
Carlo simulations we neglected this uncertainty. The second systematic uncertainty 



is due to the inability of the GEANT3-based CLAS simulation package GSIM [56 1 
to reproduce the CLAS response to electron tracks at different angles and momenta. 
To estimate this effect, we treated the six CLAS sectors as independent spectrome- 
ters. The normalized event yield measured in each sector was compared separately 
to the simulations, as shown in Fig. [71 The observed differences were compared 
sector-by-sector to remove uncertainties due to the event generator model. From 
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this comparison, we obtained an average systematic uncertainty varying for differ- 
ent data sets from 5 to 8%. 



The e^e -pair production correction introduces a systematic uncertainty due to 
the parametrization involved in the calculations. We assumed 30% uncertainty on 



the inclusive e+ cross section [47[ and 10% on the inelastic electron-carbon cross 
section [1,|4§|. The systematic uncertainty is kinematics dependent, but since the 
correction is very close to unity, the systematic uncertainty associated with it never 
exceeds 5%. 

e~-identification is based on the CC signal. The correlation of the CC to a particle 
track is dependent on the magnitude of the CC signal. Therefore, the systematic 
uncertainty is related to the mean number of photoelectrons in each kinematic 
point. The inefficiency of the CC signal is given by [1 — exp {—Nphe)], where Nphe 
is the mean number of photoelectrons. We took this probability as an estimate of 
the systematic uncertainty due to e~-identification. 

The systematic uncertainties of the radiative corrections (RC) were estimated by 
changing the various parametrizations used in the RC evaluation, and computing 
the differences with the standard parametrization. In particular, we changed the 
elastic ^^C form-factor by 10%, the quasi-elastic cross section by [10 -I- 40/(1 -I- 
Q"^ /A'P)] % and the inelastic cross section by 10%. The three obtained differences 
were combined in quadrature to extract the final systematic uncertainty of the 
radiative corrections. 

Since we applied corrections to the measured e~ momentum, we estimated the 
related systematic uncertainty. The momentum correction that brings the hy- 
drogen elastic peak to the correct position results in an average momentum shift 
of 0.003Ap/p for = 1.161 GeV and of 0.0035Ap/p for = 2.261 GeV. For 
Eo = 4.461 GcV we assumed that the possible momentum shift can be slightly 
larger than the two measured values for the lower energies, 0.004Ap/p. We shifted 
the entire data spectra by the abovementioned amounts to evaluate the relative 
change of the cross section in each bin. This modified relative change gave us 
an estimate of the uncertainty on the measured cross section. For the two lower 
beam energies, these relative shifts were suppressed by the factor accounting for 
the fixed value of the momentum at the hydrogen elastic peak position. This fac- 
tor, which reduced the uncertainty in the vicinity of the hydrogen elastic peak, 
was parametrized as [1 — exp {—{v — Q^/2Af)^/ (2(7^))], where the width parameter 
CTp was taken to be 0.3 GeV for the i?o = 1.161 GcV data and 0.4 GeV for the 
Eq = 2.261 GeV data. For the Eq — 4.461 GeV data set, the relative change in 
the cross section was taken directly as an estimate of the momentum measurement 
uncertainty. 

Thanks to recent measurements oi R — ^ performed at Jefferson Lab [s^], the 
precision of this quantity is greatly improved in our kinematic domain. Moreover, 
the structure function F2 in the interesting kinematic range is not very sensitive to 
the value of R. In fact, even a 100% systematic uncertainty on R gives only a few 
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percent uncertainty on F2. 

AF2, 1-e Ai? ^ ,AR 

Uncertainties of R given in Ref. 7] were propagated to F2 and the corresponding 
systematic uncertainties were always lower than 5%. 

Ah systematic uncertainties mentioned above and listed in Table [T] were combined in 
quadrature to give the total systematic uncertainty. When data sets from different beam 
energies overlapped, the systematic uncertainties were calculated as the central values 
using a weighted-average technique. 



Table 1: Average systematic uncertainties for the three data sets for each of the different sources. 



Source of uncertainty 


Data set 


Eo = 1.161 GeV 

[%] 


Eq ^ 2.261 GeV 

[%] 


Eo = 4.461 GeV 

[%] 


Normalization 


3.2 


1 


1 


Efficiency 


3 


3 


3 


Acceptance 


6.6 


4.7 


7.7 


e"'"e~-pair production 





0.001 


0.05 


identification 


0.04 


0.4 


0.2 


Radiative correction 


1.3 


1.7 


5 


Momentum measurement 


0.9 


0.03 


12 


Uncertainty on i? = — 


2 


1.2 


0.5 


Total 


8.5 


6.5 


20 



3. 7. Moments of the Structure Function F2 

The non-zero mass of the target leads to an undesirable mixing between QCD oper- 
ators of different spin in the OPE. To avoid this mixing in our analysis, we use Nacht- 
mann (Hlj moments M^^^((5^) instead of the usual (massless) Cornwall-Norton moments. 
The Nachtmann moments are defined as follows: 



dx^^F2{x,Q^) 



3 -I- 3(n + l)r + n{n + 2)r'^ 



(16) 



where r = -^/l + 4A'Px'^/Q^, M is the nucleon mass and ^ = 2x/ (1 + r). In the Bjorken 
limit IvP IQ'^ — > 0, the Nachtmann and Cornwall- Norton definitions coincide. 

The carbon structure function moments were evaluated according to the method 
described in Refs. @, 0|- The CLAS data were combined with the remaining world data 
on the structure function F2, along with the inclusive cross section data from Refs. 



10, 
35| 



:tion ^2iaion^__^ 

[111, [1^, 1131 M 1151, M [il H 111, iHliI, 

(see Fig.H]). The Q^.range of the CLAS data, from 0.175 to 4.95 (GeV/c)^, was divided 
into bins of width AQ^ = 0.05-0.1 (GeV/c)^. Within each bin, the world data were 
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Figure 8: The carbon structure function F2{x,Q'^) per nucleon at six different values. The open 
symbols represent experimental data obtained in the present analysis (triangles for Eq = 1.161 GeV, 
squares for Eq = 2.261 GeV and circles for Eq = 4.461 GeV) with systematic uncertainties indicated by 
the hatched areas. The full stars show data from previous experiments [l( lllll.fT3 . llSllTS, Il5l. [T^. llTtlTg. 
[Tgl . [20I . [21] . l22fl . The solid line represents the parametrization from Refs. [g. l48ll . 
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Figure 9: Integrands of the Nachtmann moments {I(M„)) at = 0.825 (GeV/c)^: circles represent the 
integrand of A/2; squares show the integrand of M4; triangles show the integrand of Ma; crosses show 
the integrand of Mg. 



shifted to the central bin value Qq, using the fit of F2{x,Q^) from Refs. [8|,|48[. The 
integrals of the data over x were performed numerically using the standard trapezoidal 
method TRAPER [H^. As an example, Fig. IH] shows the integrands of the first four 
moments as a function of x at fixed . The significance of the large x region for various 
moments can clearly be seen. 

As in Refs. 6, 7], the world data at above 5 (GeV/c)^ were analyzed in the same 
way as described above, but with a different bin size. The bin size was chosen to 
provide sufficient a;-coverage for most of the bins {AQ^/Q^ = 5 %). The results, 
together with their statistical and systematic uncertainties, are shown in Fig. [TU and 
reported in Table (2] A comparison of Fig. [TO] to the corresponding figures in Refs. 1, 01 
shows the lack of the carbon data in the interval from 5 to 40 (GeV/c)^. This 
kinematic interval is very important because here moments reach the scaling regime. 
The systematic uncertainty consists (see for example. Fig. 1 111) of experimental un- 



certainties in our data and the data given in Refs. fmjll^^ Ei, H [H, [13, H 



li, M, M, m El, m, M, M,M,M,M,M,M, M, IlIllE 



and uncertainties in the 
The first type of uncertainties was accounted for in the moment 
To estimate them we had to account for the inclusion of many 



evaluation procedure, 
evaluation point-by-point, 
data sets measured in different laboratories with different detectors. In the present anal- 
ysis, we assumed that the different experiments are independent and, therefore, only the 
systematic uncertainties within a given data set are correlated. The uncertainties in the 
evaluation procedure were obtained as described previously in Refs. @, 0|- The lowest 
moment, n — 2, includes a systematic uncertainty due to the low-x extrapolation, which 
is negligible for larger n. To estimate this uncertainty we modified the parametrization 
from Refs. [1, E^, used for the low-x extrapolation. The modification, evaluated by a 
fit to the available carbon F2 data in the low- a; region, was given by the multiplicative 
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Figure 10: The Nachtmann moments extracted from the world data, including the new CLAS results. 
Uncertainties are statistical only and since their values do not exceed a few percent (see Table [JJl, they 
are smaller than the symbol size. 
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Figure 11: Uncertainties of the Nachtmann moment AI2 in percent. The lower cross-hatched area 
represents the statistical uncertainties. The left-hatched area represents the systematic uncertainties. 
The right-hatched area represents the low- 3; extrapolation uncertainty. 



factor: 1 -I- (0.38 -I- 0MA/Q^)x°-^{1 - At = 1 (GeV/c)^ the relative change in 

the value of F2 parametrization ranges from 2% at x '--^ 10^^ up to 18% at x ~ 10^^. 
Then we compared the moments from the Table [5] with ones obtained by using the mod- 
ified low-x extrapolation and we took the difference as the estimate of the systematic 
uncertainty. As one can see from Fig. 1111 systematic uncertainties dominate at n = 2, 
while the statistical uncertainties become comparable at larger n. 

4. Phenomenological Twists Expansion 

The Q^-distributions of the n = 2,4,6 and 8 moments were fit using the method 
developed in Refs. d, IH, [13 with the following parametrization: 

Mnias, g') = Vn{as) + HT,,{as, Q^) , (17) 

where for n > 4 the leading-twist term contained only one unknown parameter due to the 
dominant non-singlet contribution, while for n — 2, both singlet and non-singlet terms 
were considered (for details see Refs. 0, 0])- 

The higher- twist contribution is given by [6l| : 

adding 2A^ free parameters AnT and 7m- , with N being the number of higher-twist terms 
considered. Previous analyses [1, 0, [m] showed that it was necessary to take N > 2 
when n > 4, while for n = 2 it was sufficient to have N > 1. In principle, the number 
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(18) 



Tabic 2: The Nachtmann moments for n = 2,4,6 and 8 evaluated in the interval 0.15 < Q"^ < 150 (GeV/c)'^. The moments are labeled with an 
asterisk when the contribution to the integral by the experimental data is between 50% and 70%. All the others were evaluated with more than 70% 
data coverage. The data are reported together with the statistical and systematic uncertainties, the third uncertainty for n = 2 is due to the low- a; 
extrapolation. 



Q K^ev/cj J 


M2{Q jxiU 


M4[Q jxiU 


Me[CJ jxiU 


i\/f /"/^^ r\— 
Ms{Q JxiO 


0.175 


1.602 ± 0.003 ± 0.126 ± 0.004 


1.707 ± 0.004 ± 0.133 


0.222 ± 0.001 ± 0.018 


0.302 ± 0.001 ± 0.027 


0.225 


1.915 ± 0.003 ± 0.223 ± 0.003 


2.428 ± 0.007 ± 0.179 


0.395 ± 0.002 ± 0.029 


0.687 ± 0.003 ± 0.051 


0.275 


1.923 ± 0.003 ± 0.645 ± 0.001 


2.715 ± 0.005 ± 0.303 


0.500 ± 0.001 ± 0.045 


0.990 ± 0.002 ± 0.092 


0.325 


1.991 ± 0.004 ± 0.256 ± 0.003 


3.063 ± 0.006 ± 0.198 


0.635 ± 0.001 ± 0.040 


1.424 ± 0.004 ± 0.100 


0.375 


1.985 ± 0.003 ± 0.592 ± 0.005 


3.273 ± 0.007 ± 0.271 


0.746 ± 0.002 ± 0.044 


1.854 ± 0.006 ± 0.111 


0.425 




3.452 ± 0.008 ± 0.215 


0.847 ± 0.002 ± 0.057 


2.298 ± 0.008 ± 0.180 


0.475 


2.067 ± 0.009 ± 0.145 ± 0.000 


3.544 ± 0.009 ± 0.219 


0.926 ± 0.003 ± 0.059 


2.707 ± 0.011 ± 0.188 


0.525 


2.091 ± 0.003 ± 0.391 ± 0.005 


3.608 ± 0.008 ± 0.231 


0.990 ± 0.003 ± 0.062 


3.073 ± 0.011 ± 0.222 


0.575 


2.090 ± 0.002 ± 0.176 ± 0.005 


3.667 ± 0.008 ± 0.194 


1.051 ± 0.003 ± 0.057 


3.435 ± 0.013 ± 0.216 


0.625 


2.061 ± 0.002 ± 0.150 ± 0.002 


3.636 ± 0.009 ± 0.209 


1.085 ± 0.004 ± 0.065 


3.730 ± 0.017 ± 0.260 


0.675 


2.011 ± 0.003 ± 0.254 ± 0.002 


3.577 ± 0.010 ± 0.225 


1.102 ± 0.005 ± 0.067 


3.949 ± 0.024 ± 0.276 


0.725 


1.977 ± 0.002 ± 0.326 ± 0.001 


3.537 ± 0.010 ± 0.242 


1.125 ± 0.005 ± 0.072 


4.208 ± 0.028 ± 0.313 


0.775 




3.488 ± 0.011 ± 0.213 


1.139 ± 0.006 ± 0.069 


4.417 ± 0.034 ± 0.324 


0.825 




3.418 ± 0.009 ± 0.205 


1.139 ± 0.004 ± 0.072 


4.545 ± 0.024 ± 0.332 


0.875 




3.357 ± 0.008 ± 0.201 


1.148 ± 0.005 ± 0.074 


4.741 ± 0.027 ± 0.369 


0.925 




3.295 ± 0.008 ± 0.327 


1.143 ± 0.005 ± 0.082 


4.835 ± 0.027 ± 0.405 


0.975 




3.230 ± 0.009 ± 0.229 


1.144 ± 0.005 ± 0.066 


4.976 ± 0.030 ± 0.289 



i(GcV/c)-J 


Maig-^jxlO-^ 


M4(O^)xl0-^ 


M6(Q")xlO-^ 


Ms{Q')xlO-^ 


1.050 




3.097 ± 0.006 ± 0.205 


1.109 ± 0.004 ± 0.065 


4.936 ± 0.023 ± 0.277 


1.150 




2.911 ± 0.007 ± 0.234 


1.058 ± 0.004 ± 0.077 


4.911 ± 0.029 ± 0.402 


1.250 




2.799 ± 0.007 ± 0.259 


1.035 ± 0.005 ± 0.081 


4.957 ± 0.035 ± 0.429 


1.350 


1.672 ± 0.010 ± 0.130 ± 0.004 


2.673 ± 0.007 ± 0.221 


0.984 ± 0.004 ± 0.080 


4.778 ± 0.034 ± 0.421 


1.450 


*1.622 ± 0.002 ± 0.141 ± 0.017 


2.541 ± 0.008 ± 0.215 


0.938 ± 0.005 ± 0.076 


4.627 ± 0.034 ± 0.385 


1.550 




2.443 ± 0.008 ± 0.254 


0.889 ± 0.005 ± 0.081 


4.408 ± 0.039 ± 0.416 


1.650 




2.343 ± 0.008 ± 0.267 


0.847 ± 0.005 ± 0.083 


4.232 ± 0.036 ± 0.437 


1.750 




2.250 ± 0.009 ± 0.238 


0.805 ± 0.005 ± 0.089 


4.033 ± 0.038 ± 0.489 


1.850 




2.174 ± 0.010 ± 0.244 


0.766 ± 0.005 ± 0.096 


3.844 ± 0.040 ± 0.575 


1.950 


*1.540 ± 0.006 ± 0.114 ± 0.014 


2.067 ± 0.010 ± 0.221 


0.718 ± 0.005 ± 0.085 


3.582 ± 0.036 ± 0.473 


2.050 




2.068 ± 0.012 ± 0.235 


0.725 ± 0.007 ± 0.096 


3.731 ± 0.047 ± 0.570 


2.150 




2.027 ± 0.012 ± 0.257 


0.696 ± 0.007 ± 0.098 


3.533 ± 0.050 ± 0.586 


2.250 




1.974 ± 0.013 ± 0.228 


0.687 ± 0.008 ± 0.105 


3.560 ± 0.066 ± 0.663 


2.350 




1.891 ± 0.013 ± 0.225 


0.642 ± 0.009 ± 0.113 


3.300 ± 0.080 ± 0.791 


2.450 




1.901 ± 0.013 ± 0.424 


0.630 ± 0.009 ± 0.130 


3.180 ± 0.087 ± 0.814 


2.550 






0.605 ± 0.008 ± 0.096 


3.042 ± 0.065 ± 0.614 


2.650 






0.590 ± 0.009 ± 0.104 


2.983 ± 0.082 ± 0.721 


2.750 






0.585 ± 0.009 ± 0.092 


2.986 ± 0.075 ± 0.615 


2.850 






0.557 ± 0.008 ± 0.104 


2.741 ± 0.069 ± 0.643 


2.950 






0.557 ± 0.009 ± 0.094 


2.831 ± 0.082 ± 0.603 



i(GcV/c)^J 


M2(Q^)xl0-^ 


M4(Q^)xlO-^ 


A^6(Q')xlO-^ 


Ms{Q'')xlO-'' 


3.050 






0.540 ± 0.009 ± 0.095 


2.717 ± 0.075 ± 0.611 


3.150 






0.520 ± 0.008 ± 0.092 


2.616 ± 0.075 ± 0.601 


3.250 






0.513 ± 0.008 ± 0.103 


2.536 ± 0.070 ± 0.642 


3.350 








2.597 ± 0.068 ± 0.655 


3.450 








2.409 ± 0.069 ± 0.579 


3.550 








2.437 ± 0.056 ± 0.747 


3.650 








2.295 ± 0.078 ± 0.765 


3.750 








2.354 ± 0.072 ± 0.613 


3.850 








2.137 ± 0.063 ± 0.589 


3.950 








2.170 ± 0.064 ± 0.599 


4.050 








2.132 ± 0.072 ± 0.557 


4.150 








2.038 ± 0.071 ± 0.692 


5.125 


*1.457 ± 0.004 ± 0.047 ± 0.007 








62.500 




0.920 ± 0.007 ± 0.023 


0.189 ± 0.002 ± 0.002 


0.608 ± 0.009 ± 0.009 


83.250 




0.877 ± 0.009 ± 0.016 


0.177 ± 0.003 ± 0.001 


0.561 ± 0.010 ± 0.007 


87.000 




*0.874 ± 0.007 ± 0.026 


0.176 ± 0.002 ± 0.007 




96.000 




*0.824 ± 0.011 ± 0.030 






151.500 




*0.821 ± 0.015 ± 0.018 


0.162 ± 0.005 ± 0.001 


0.512 ± 0.020 ± 0.003 



Table 3: Extracted parameters of the twist expansion at the reference scale /x^ = 10 (GeV/c)^. The first 
uncertainty has a statistical origin and is obtained from a MINOS [sol minimization procedure, while 
the second uncertainty is the systematic one described in text. The contribution of twist-6 to M2 was 
too small to be extracted by the present procedure. 



n 


2 


4 


6 


8 




0.1411±0.0005±0.0050 


(1.206±0.008±0.051) x 10"^ 


(2.79±0.05±0.07)xl0"^ 


(9.7±0.5±0.3)xl0"*' 


A-n4 


(0.6±1.8±l.l)xl0"'* 


(2.14±0.02±2.68)xl0"'' 


(2.22±0.03±3.25)xl0"'* 


(2.9±0.1±2.7)xl0"'' 


7n4 


8±3±10 


4.95±0.01±1.65 


5.21±0.02±4.75 


3.7±0.1±2.1 


A-n6 




(-4.65±0.04±5.07)xl0"'^ 


(-6.2±0.1±5.5)xl0"^ 


(-7.5±0.4±4.1)xl0"''* 


7n6 




3.83±0.01±1.58 


3.84±0.02±2.21 


2.5±0.1±1.9 



of higher- twist terms N should be infinite, but as was shown in Ref. [6^, aheady the 
minimum number of terms (i.e. N = 1 for n = 2 and N = 2 when n > 4) completely 
described the total higher-twist contribution. Adding further terms does not change the 
obtained separation between leading and higher twists. Hence, in the present analysis 
we used = 1 for n = 2 and N = 2 when n > 4. 

The higher-twist parameters 7„t- are effectively the LO anomalous dimensions of the 
perturbative Wilson coefficients Em-i^J■r, ^J■f i cts) in Eq. [21 while the parameters Anr are 
proportional to the matrix elements of QCD operators OnrifJ-)- Higher- twist operators 
get mixed and only leading-twist parameters An2 can be directly related to the operators 
0„2(/^) for r > 2 since pQCD calculations of "fnr are at present unavailable. 

Based on the above, the n-th moment (see Ea.ll8|) for n > A has five unknown parame- 
ters: the leading-twist parameter A„2 and the higher-twist parameters A„4, 7„4, A„g and 
7„6- These parameters were simultaneously determined from a x^-minimization proce- 
dure (MINUIT [63]) over the allowed range. The statistical uncertainties of the exper- 
imental moments were used by MINOS [6^ to obtain statistical uncertainties on the ex- 
tracted parameters. Their systematic uncertainties were obtained by adding/subtracting 
the systematic uncertainties to the experimental moments and by repeating the twist 
extraction procedure. 

In the case of n = 2, the sum of the non-singlet and singlet terms at the leading twist 
was considered, therefore adding the leading-twist parameter due to the giuon moment 
at the reference scale. Moreover, due to the vanishing contribution of the higher twists 
in M2 (see Fig. [TS]). we considered the twist-4 term only. 

The fit results are shown in Fig.[T21 while in Table [31 we report the parameter values 
obtained at the reference scale fj,^ = 10 (GeV/c)^. In addition, the extracted leading- 
twist contribution is reported in Table [H 

5. Comparison to Deuteron 

Once the separation between the leading and higher-twist terms is obtained, we can 
compare them to the corresponding terms in the deuteron. The deuteron, being a loosely 
bound system, is considered here as representing structure function moments of an almost 
free proton and neutron. 

In Fig. [T2] the ratio of the leading-twist moments of carbon-to-deuteron is shown 
with its statistical and systematic uncertainties. It is compared to the ratio of the 
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Figure 12: Results of the twist analysis. The squares represent the Nachtmann moments obtained in 
this analysis. The solid line is the fit to the moments using Eg. 1 171 with the parameters listed in Table [51 
The leading-twist (LT), twist-4 (TW-4), twist-6 (TW-6) and the total higher twist (HT) contributions 
to the fit are indicated. 



26 



Table 4: The extracted leading- twist contribution r]n{Q^) shown in Fig. 1121 reported with systematic 
uncertainties. 



Q' [(GeV/c)^] 


772 (Q')xlG-^ 


774(Q^)xlO-" 


7?6(Q")xlO-^ 


778(g^)xlO-^ 


1.025 


1.694 ± 0.004 


2.54 ± 0.02 


10.88 ± 0.20 


7.97 ± 0.39 


1.075 


1.685 ± 0.004 


2.47 ± 0.02 


10.03 ± 0.18 


6.72 ± 0.33 


1.125 


1.675 ± 0.003 


2.40 ± 0.02 


9.35 ± 0.17 


5.85 ± 0.28 


1.175 


1.666 ± 0.003 


2.34 ± 0.02 


8.80 ± 0.16 


5.21 ± 0.25 


1.225 


1.657 ± 0.003 


2.29 ± 0.02 


8.33 ± 0.15 


4.71 ± 0.23 


1.275 


1.648 ± 0.003 


2.24 ± 0.02 


7.93 ± 0.14 


4.32 ± 0.21 


1.325 


1.640 ± 0.003 


2.19 ± 0.02 


7.59 ± 0.14 


4.00 ± 0.19 


1.375 


1.632 ± 0.003 


2.15 ± 0.01 


7.29 ± 0.13 


3.74 ± 0.18 


1.425 


1.624 ± 0.002 


2.11 ± 0.01 


7.03 ± 0.13 


3.52 ± 0.17 


1.475 


1.616 ± 0.002 


2.08 ± 0.01 


6.80 ± 0.12 


3.33 ± 0.16 


1.525 


1.609 ± 0.002 


2.05 ± 0.01 


6.59 ± 0.12 


3.17 ± 0.15 


1.575 


1.602 ± 0.002 


2.02 ± 0.01 


6.40 ± 0.11 


3.02 ± 0.14 


1.625 


1.596 ± 0.002 


1.99 ± 0.01 


6.23 ± 0.11 


2.90 ± 0.14 


1.675 


1.590 ± 0.002 


1.96 ± 0.01 


6.07 ± 0.11 


2.79 ± 0.13 


1.725 


1.584 ± 0.002 


1.94 ± 0.01 


5.93 ± 0.11 


2.69 ± 0.13 


1.775 


1.578 ± 0.002 


1.91 ± 0.01 


5.80 ± 0.10 


2.60 ± 0.12 


1.825 


1.572 ± 0.002 


1.89 ± 0.01 


5.68 ± 0.10 


2.51 ± 0.12 


1.875 


1.567 ± 0.002 


1.87 ± 0.01 


5.56 ± 0.10 


2.44 ± 0.12 


1.925 


1.562 ± 0.002 


1.85 ± 0.01 


5.46 ± 0.10 


2.37 ± 0.11 


1.975 


1.557 ± 0.002 


1.83 ± 0.01 


5.36 ± 0.09 


2.31 ± 0.11 


2.025 


1.552 ± 0.002 


1.81 ± 0.01 


5.27 ± 0.09 


2.25 ± 0.11 


2.075 


1.548 ± 0.002 


1.80 ± 0.01 


5.18 ± 0.09 


2.20 ± 0.10 


2.125 


1.543 ± 0.002 


1.78 ± 0.01 


5.10 ± 0.09 


2.15 ± 0.10 


2.175 


1.539 ± 0.002 


1.76 ± 0.01 


5.02 ± 0.09 


2.10 ± 0.10 


2.225 


1.535 ± 0.002 


1.75 ± 0.01 


4.95 ± 0.09 


2.06 ± 0.10 


2.275 


1.532 ± 0.002 


1.74 ± 0.01 


4.89 ± 0.09 


2.02 ± 0.10 


2.325 


1.529 ± 0.002 


1.72 ± 0.01 


4.83 ± 0.08 


1.99 ± 0.09 


2.375 


1.527 ± 0.002 


1.71 ± 0.01 


4.78 ± 0.08 


1.96 ± 0.09 
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Q [(GeV/c)^J 


ri2{Q )xlO 


t' /~\'Z\ 1 r\ — 'Z 


ryelQ jxlO 


r/8((y jxlO 


2.425 


1.524 ± 0.002 


1.70 ± 0.01 


A TO 1 A AO 

4.73 ± 0.08 


1 AO 1 A AA 

1.93 ± 0.09 


2.475 


1.522 ± 0.002 


1 /^A 1 A Al 

1.69 ± 0.01 


A /^O 1 A AO 

4.68 ± 0.08 


1 AA 1 A AA 

1.90 ± 0.09 


2.525 


1.520 ± 0.002 


1 O 1 A A 1 

1.68 ± 0.01 


A i~t 1 A A O 

4.63 ± 0.08 


1 OO 1 A A A 

1.88 ± 0.09 


2.575 


1.518 ± 0.002 


1 T 1 A A 1 

1.67 ± 0.01 


A r n 1 A AO 

4.59 ± 0.08 


1 O r" 1 A A A 

1.85 ± 0.09 


2.625 


1.516 ± 0.002 


1 /^/^ 1 A A"! 

1.66 ± 0.01 


A f A 1 A AO 

4.54 ± 0.08 


1 OO 1 A AA 

1.83 ± 0.09 


2.675 


1.514 ± 0.002 


1 r" 1 A A 1 

1.65 ± 0.01 


A r f\ 1 A A O 

4.50 ± 0.08 


1 O A 1 A A O 

1.80 ± 0.08 


2.725 


1.512 ± 0.002 


1/^/1 1 A A 1 

1.64 ± 0.01 


A A 1 A AO 

4.46 ± 0.08 


1 'TO 1 A AO 

1.78 ± 0.08 


2.775 


1.510 ± 0.002 


1/^0 1 A A 1 

1.63 ± 0.01 


A A C\ 1 A AO 

4.42 ± 0.08 


1 1 A AO 

1.76 ± 0.08 


2.825 


1.508 ± 0.002 


1 A 1 A A 1 

1.62 ± 0.01 


< O A 1 A A O 

4.39 ± 0.08 


'A 'T A 1 A A O 

1.74 ± 0.08 


2.875 


1.506 ± 0.002 


1/^1 1 A A 1 

1.61 ± 0.01 


^ O P" 1 A AO 

4.35 ± 0.08 


1 TO 1 A AO 

1.72 ± 0.08 


o no c 


1.504 ± 0.003 


1.61 ± 0.01 


4.62 ± O.Oo 


1.70 ± O.Oo 


2.975 


1.503 ± 0.003 


1.60 ± 0.01 


4.28 ± 0.07 


1.69 ± 0.08 


3.025 


1.501 ± 0.003 


1.59 ± 0.01 


4.25 ± 0.07 


1.67 ± 0.08 


3.075 


1.499 ± 0.003 


1.58 ± 0.01 


4.22 ± 0.07 


1.65 ± 0.08 


3.125 


1.498 ± 0.003 


1.58 ± 0.01 


4.19 ± 0.07 


1.64 ± 0.08 


3.175 


1.496 ± 0.003 


1.57 ± 0.01 


4.16 ± 0.07 


1.62 ± 0.08 


3.225 


1.495 ± 0.003 


1.56 ± 0.01 


4.13 ± 0.07 


1.61 ± 0.07 


3.275 


1.493 ± 0.003 


1.56 ± 0.01 


4.10 ± 0.07 


1.59 ± 0.07 


3.325 


1.492 ± 0.003 


1.55 ± 0.01 


4.08 ± 0.07 


1.58 ± 0.07 


3.375 


1.490 ± 0.003 


1.54 ± 0.01 


4.05 ± 0.07 


1.56 ± 0.07 


3.425 


1.489 ± 0.003 


1.54 ± 0.01 


4.02 ± 0.07 


1.55 ± 0.07 


3.475 


1.488 ± 0.003 


1.53 ± 0.01 


4.00 ± 0.07 


1.54 ± 0.07 


3.525 


1.486 ± 0.003 


1.52 ± 0.01 


3.98 ± 0.07 


1.53 ± 0.07 


3.575 


1.485 ± 0.003 


1.52 ± 0.01 


3.95 ± 0.07 


1.51 ± 0.07 


3.625 


1.484 ± 0.003 


1.51 ± 0.01 


3.93 ± 0.07 


1.50 ± 0.07 


3.675 


1.482 ± 0.003 


1.51 ± 0.01 


3.91 ± 0.07 


1.49 ± 0.07 


3.725 


1.481 ± 0.003 


1.50 ± 0.01 


3.89 ± 0.07 


1.48 ± 0.07 


3.775 


1.480 ± 0.003 


1.50 ± 0.01 


3.87 ± 0.07 


1.47 ± 0.07 


3.825 


1.479 ± 0.003 


1.49 ± 0.01 


3.84 ± 0.07 


1.46 ± 0.07 


3.875 


1.478 ± 0.003 


1.49 ± 0.01 


3.82 ± 0.07 


1.45 ± 0.07 
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Q [(GeV/c)^J 


ri2{Q )xlO 


t' /~\'Z\ 1 r\ — 'Z 


ryelQ jxlO 


r/8((y jxlO 


3.925 


1.477 ± 0.003 


1 AC) 1 /"» A 1 

1.48 ± 0.01 


O O 1 1 A A'? 

3.81 ± 0.07 


-1 A A 1 A AT 

1.44 ± 0.07 


3.975 


1.476 ± 0.003 


1.48 ± 0.01 


O ^A 1 A A^ 

3.79 ± 0.07 


1 /I O 1 A AT 

1.43 ± 0.07 


4.025 


1.474 ± 0.003 


1.47 ± 0.01 


3.77 ± 0.07 


1 /I O 1 A AT 

1.42 ± 0.07 


4.075 


1.473 ± 0.003 


1.47 ± 0.01 


3.75 ± 0.06 


1/11 1 A AT 

1.41 ± 0.07 


4.125 


1.472 ± 0.003 


1.46 ± 0.01 


O 1 A A/^ 

3.73 ± 0.06 


-t Af\ 1 A A/^ 

1.40 ± 0.06 


4.175 


1.471 ± 0.003 


1 A 1 Ci A 1 

1.46 ± 0.01 


O T1 1 A Ar* 

3.71 ± 0.06 


-^ A r\ 1 A A/^ 

1.40 ± 0.06 


4.225 


1.470 ± 0.003 


1.45 ± 0.01 


O ^A 1 A A/^ 

3.70 ± 0.06 


1 OA 1 A A/^ 

1.39 ± 0.06 


4.275 


1.469 ± 0.004 


1 A f 1 A A 1 

1.45 ± 0.01 


rj /^o 1 A A/' 

3.68 ± 0.06 


1 OO 1 A A/^ 

1.38 ± 0.06 


4.325 


1 A C> 1 r\ r\f\ A 

1.468 ± 0.004 


1 /IP" 1 A A 1 

1.45 ± 0.01 


^ O O 1 A A*^ 

3.66 ± 0.06 


1 OT 1 A A/^ 

1.37 ± 0.06 


4.375 


1.467 ± 0.004 


1 A A 1 A A 1 

1.44 ± 0.01 


o r 1 A A/^ 

3.65 ± 0.06 


1 O/^ 1 A A/^ 

1.36 ± 0.06 


A /IOC 

4.4Z0 


1.467 ± 0.004 


1.44 ± 0.01 


3.63 ± 0.06 


l.oo ± 0.06 


4.475 


1.466 ± 0.004 


1.43 ± 0.01 


3.62 ± 0.06 


1.35 ± 0.06 


4.525 


1.465 ± 0.004 


1.43 ± 0.01 


3.60 ± 0.06 


1.34 ± 0.06 


4.575 


1.464 ± 0.004 


1.43 ± 0.01 


3.59 ± 0.06 


1.33 ± 0.06 


4.625 


1.463 ± 0.004 


1.42 ± 0.01 


3.57 ± 0.06 


1.33 ± 0.06 


4.675 


1.462 ± 0.004 


1.42 ± 0.01 


3.56 ± 0.06 


1.32 ± 0.06 


4.725 


1.461 ± 0.004 


1.41 ± 0.01 


3.54 ± 0.06 


1.31 ± 0.06 


4.775 


1.460 ± 0.004 


1.41 ± 0.01 


3.53 ± 0.06 


1.31 ± 0.06 


4.825 


1.460 ± 0.004 


1.41 ± 0.01 


3.52 ± 0.06 


1.30 ± 0.06 


4.875 


1.459 ± 0.004 


1.40 ± 0.01 


3.50 ± 0.06 


1.29 ± 0.06 


4.925 


1.458 ± 0.004 


1.40 ± 0.01 


3.49 ± 0.06 


1.29 ± 0.06 


4.975 


1.457 ± 0.004 


1.40 ± 0.01 


3.48 ± 0.06 


1.28 ± 0.06 


5.025 


1.456 ± 0.004 


1.39 ± 0.01 


3.47 ± 0.06 


1.27 ± 0.06 


5.075 


1.456 ± 0.004 


1.39 ± 0.01 


3.45 ± 0.06 


1.27 ± 0.06 


5.125 


1.455 ± 0.004 


1.39 ± 0.01 


3.44 ± 0.06 


1.26 ± 0.06 


5.275 


1.453 ± 0.004 


1.38 ± 0.01 


3.41 ± 0.06 


1.25 ± 0.06 


5.325 


1.452 ± 0.004 


1.37 ± 0.01 


3.40 ± 0.06 


1.24 ± 0.06 


5.375 


1.451 ± 0.004 


1.37 ± 0.01 


3.38 ± 0.06 


1.24 ± 0.06 


5.475 


1.450 ± 0.004 


1.37 ± 0.01 


3.36 ± 0.06 


1.23 ± 0.06 


5.525 


1.449 ± 0.004 


1.36 ± 0.01 


3.35 ± 0.06 


1.22 ± 0.06 
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[{Ge\/cf] 


r/2(0")xlO-^ 


r;4(Q")xl0-^ 


%(Q")xl0-^ 


r?8(Q")xlO-^ 


5.625 


1.448 ± 0.004 


1.36 ± 0.009 


3.33 ± 0.06 


1.21 ± 0.06 


5.675 


1.447 ± 0.004 


1.35 ± 0.009 


3.32 ± 0.06 


1.21 ± 0.05 


5.725 


1.447 ± 0.004 


1.35 ± 0.009 


3.31 ± 0.06 


1.20 ± 0.05 


5.955 


1.444 ± 0.004 


1.34 ± 0.009 


3.27 ± 0.06 


1.18 ± 0.05 


6.915 


1.434 ± 0.005 


1.30 ± 0.009 


3.11 ± 0.05 


1.11 ± 0.05 


7.267 


1.431 ± 0.005 


1.28 ± 0.008 


3.06 ± 0.05 


1.09 ± 0.05 


7.630 


1.427 ± 0.005 


1.27 ± 0.008 


3.02 ± 0.05 


1.07 ± 0.05 


8.021 


1.424 ± 0.005 


1.26 ± 0.008 


2.97 ± 0.05 


1.05 ± 0.05 


8.847 


1.418 ± 0.005 


1.23 ± 0.008 


2.89 ± 0.05 


1.01 ± 0.05 


9.775 


1.413 ± 0.005 


1.21 ± 0.008 


2.80 ± 0.05 


0.97 ± 0.04 


10.267 


1.410 ± 0.005 


1.20 ± 0.008 


2.77 ± 0.05 


0.95 ± 0.04 


10.762 


1.407 ± 0.005 


1.18 ± 0.008 


2.73 ± 0.05 


0.94 ± 0.04 


11.344 


1.405 ± 0.005 


1.17 ± 0.008 


2.69 ± 0.05 


0.92 ± 0.04 


12.580 


1.399 ± 0.006 


1.15 ± 0.008 


2.62 ± 0.04 


0.89 ± 0.04 


13.238 


1.397 ± 0.006 


1.14 ± 0.007 


2.58 ± 0.04 


0.88 ± 0.04 


14.689 


1.392 ± 0.006 


1.12 ± 0.007 


2.52 ± 0.04 


0.85 ± 0.04 


17.108 


1.385 ± 0.006 


1.09 ± 0.007 


2.42 ± 0.04 


0.81 ± 0.04 


19.072 


1.380 ± 0.006 


1.07 ± 0.007 


2.36 ± 0.04 


0.78 ± 0.03 


20.108 


1.378 ± 0.006 


1.06 ± 0.007 


2.33 ± 0.04 


0.77 ± 0.03 


21.097 


1.376 ± 0.006 


1.05 ± 0.007 


2.31 ± 0.04 


0.76 ± 0.03 


24.259 


1.370 ± 0.006 


1.03 ± 0.007 


2.23 ± 0.04 


0.73 ± 0.03 


26.680 


1.369 ± 0.007 


1.01 ± 0.007 


2.19 ± 0.04 


0.71 ± 0.03 


32.500 


1.367 ± 0.007 


0.99 ± 0.006 


2.10 ± 0.03 


0.68 ± 0.03 


34.932 


1.367 ± 0.007 


0.98 ± 0.006 


2.07 ± 0.03 


0.67 ± 0.03 


36.750 


1.366 ± 0.007 


0.97 ± 0.006 


2.05 ± 0.03 


0.66 ± 0.03 


43.970 


1.365 ± 0.007 


0.94 ± 0.006 


1.99 ± 0.03 


0.63 ± 0.03 


47.440 


1.365 ± 0.007 


0.93 ± 0.006 


1.96 ± 0.03 


0.62 ± 0.03 


64.270 


1.363 ± 0.007 


0.90 ± 0.006 


1.85 ± 0.03 


0.58 ± 0.02 


75.000 


1.363 ± 0.007 


0.88 ± 0.006 


1.80 ± 0.03 


0.56 ± 0.02 


86.000 


1.362 ± 0.007 


0.87 ± 0.006 


1.76 ± 0.03 


0.55 ± 0.02 


97.690 


1.362 ± 0.007 


0.85 ± 0.005 


1.72 ± 0.03 


0.53 ± 0.02 
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corresponding structure functions obtained in Ref. 35| and satisfied the cut > 25 
(GeV/c)^, ensuring leading-twist dominance. 

Considering Fermi motion effects only, and applying the Impulse Approximation 

(lA) [ii,!!!: 

F,''ix,Q')^J dzf^{z)F^[-^,Q^) , (19) 

where f^{z) is the non-relativistic nucleon (light-cone) momentum distribution in the 
nucleus and the superscripts A and N are introduced to distinguish the nuclear (A) and 
nucleon (iV) distributions. Then, if we define the n-th moment of the distribution f^{z) 
by: 

pMa/M 

= / dzz-'^!\z) , (20) 



in the moment space, the convolution becomes the product of the moments [66[: 

M,f(Q2)= j-„^M„^(q2) (21) 

This leads to a simple relation between the carbon and deuteron moments in the lA: 



(22) 



where the nucleon structure function moments are canceled. Therefore, in the lA the ratio 
of nuclear structure function moments is reduced to the ratio of the nucleon momentum 
distribution moments. This ratio, obtained with nuclear wave functions from Refs. |48l . 



67[, is shown in Fig.[T51bv the solid line. The deviation of the data points from the curve 
at large n is the consequence of the EMC effect. Indeed, the ratio of the moments is in 
good agreement with the ratio of the structure functions taken at x values corresponding 
to the average < a; >„ for a given n defined as in Ref. (68l |. Hence, in the moment space, 
we confirm the EMC effect discovered previously in a;-space. 

A similar ratio of the carbon to deuteron moments, but for the total higher-twist 
contribution taken at = 2 (GcV/c)^, is shown in Fig. [TH The curve representing 
the Fermi motion expectation is the same as for the leading-twist case. The higher-twist 
ratio is very different from the ratio of leading twists. Surprisingly, for small n < 8, the 
total higher-twist contribution in complex nuclei is smaller than in the deuteron. Despite 
the large systematic uncertainties, the ratio rises with n almost linearly, surpassing unity 
at n « 7. This behavior is clearly not related to the Fermi motion of nucleons in the 
nucleus. Naively one may expect an overall larger total higher-twist contribution in 
mplex nuclei because of additional processes, e.g. nuclear FSI and SRC [H [H [11], 



cor 

which are less important in the deuteron. A hint of this suppression of the higher twists 
in nuclei was already observed in Ref. y, where the authors tried to connect it to the 
suppression of the nucleon excited state form- factors in nuclei (toI . Q ■ Indeed, looking at 
the resonance peaks in the free proton, deuteron and carbon, shown in Fig. 1151 one may 
conclude that they disappear in nuclei. However, the suppression factors extracted in 
Ref. would lead to a much larger difference between the carbon and deuteron moments. 
This is because the total resonance contribution to the structure function moments in the 
few (GeV/c)^ domain is very important [4]. Thus, if resonance form-factors are really 
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Figure 13: The carbon-to-deuteron ratio of the leading-twist moments (full triangles) compared with the 
ratio of corresponding structure functions obtained in Ref. 35] (open squares) taken at average <^ x 
values (the cut > 25 (GeV/c)^ is applied to avoid higher-twist contamination in the structure 
function ratio). The curve shows the theoretical expectation due to Fermi motion effects. Internal error 
bars are statistical only and the total error bars arc the statistical and systematic uncertainties summed 
in quadrature. 
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suppressed in nuclei, this has to be compensated by an almost equivalent increase of the 
non-resonant contribution to the first four moments. 

One explanation of the EMC effect as due to partial quark deconfinement for a nu- 
cleon embedded in nuclear matter was suggested within the rescaling model proposed 
in Ref. 71]. In Ref. [8|] it was also connected to the higher-twist suppression in nuclei. 
In fact, if one considers the normal QCD potential, partons inside the nucleon will in- 
teract more strongly when approaching the nucleon radius, while partons deep inside 
the nucleon will behave as free particles. This intuitive picture is the basis for the MIT 



Bag model (T^- The discussion below is a simple-minded speculation suggested by this 
picture. 

The interaction of the struck parton with the rest of the nucleon will result in the 
appearance of higher twists, but as long as the parton is located in the central region of 
the nucleon, it will propagate as a free particle, resulting in the leading-twist dominance. 
Therefore, for a given x (or alternatively n in the moment space), the value at which 
the higher-twist contribution becomes significant '^^^ considered as the inverse 

of the mean free path of the struck parton: 



/ Af2 

where the nucleon charge radius r^r and nucleon mass M define the characteristic scale. 
An increase of the nucleon radius as in the model of Ref. [tiI would lead to a shift of 
the higher-twist contribution to lower for all n. Both observed ratios of the leading 
and higher twists are compatible with such a naive picture. One can combine these 
observations together to visualize the spatial distribution of the free partons inside the 
nucleon as shown in Fig. [121 Here we defined Q^^rp used in Eq. [53] as the value at 
which the total higher-twist contribution in the given nth moment reaches 5% of the 
leading twist. As one can see in the nucleon bound inside the carbon nucleus, partons 
have a larger free propagation range, but the amount of highly energetic free partons is 
decreased, keeping the integrated strength almost constant. 



6. Conclusions 

A measurement of inclusive electron scattering by a carbon nucleus was performed in 
a wide two-dimensional range of x and Q^. This measurement improves the kinematic 
coverage of previously available data in the low to medium domain and when combined 
with the existing world data set allows a calculation of the Nachtmann moments for 
n = 2, 4, 6 and 8 over the Q^-range 0.2 — 150 (GeV/c)^. The carbon-to-deuteron ratio of 
the leading-twist F2 moments exhibits the well known EMC effect by deviating from the 
Fermi motion expectation. The deviation is compatible with that observed previously in 
X space in DIS. For the first time we obtained the carbon-to-deuteron ratio of the total 
higher-twist contributions to the F2 moments of the nucleon. Despite large systematic 
uncertainties, this ratio has a surprising behavior, increasing almost linearly with the 
moment order n, and lies well below unity for n < 7. This suppression of the higher 
twists in the nucleon bound in the nuclear matter cannot be described by the strong 
damping of the nucleon excitations in nuclei suggested in Ref. 0. The comparison 
between the relative contributions of higher twists in carbon and deuteron suggests a 
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Figure 14: The carbon-to-deuteron ratio of the higher-twists moments at = 2 (GeV/c)'^. The curve 
is the same as in Fig. 1131 The systematic uncertainties are shown by the hatched histogram. 
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Figure 15: The nucleon structure function F2 in the resonance region for proton (full triangles), 
deuteron Q (open circles) and carbon (full squares) at = 1 (GeV/c)^ (left) and = 3 (GeV/c)^ 
(right). Errors are statistical only. 
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Figure 16: Mean momentum of partons in the nucleon in the deuteron (open squares) and in carbon 
(full triangles) as a function of their mean free path defined in the text (see Eq. I23II . The curves show 
phenomenological Gaussian parametrizations for the deuteron (solid) and carbon (dashed). Errors are 
statistical only. 



wider distribution of free partons in the bound nucleon. We speculate therefore that 
the spatial shape of the nucleon, thought as a bag of free partons, is enhanced when it 
is immersed in the nuclear matter, but in contrast to the rescaling model, the overall 
probability to find a parton is conserved. 

Our analysis indicates a need for further measurements encompassing the kinematic 
region of from 5 to 40 (GeV/c)^ and large x, and > 10 (GcV/c)^ and low x. The 
first region can be explored with 12-GeV Jefferson Lab upgrade [73|, while the second 
domain would be accessible with the construction of an electron- ion collider [t^] ■ 
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